Collinear versus non-collinear magnetic order in Pd atomic clusters 



O 

o 

C 

o 

(N 



CO 



c3 



c 

o 
o 



in 

(N 
O 

o 



X 



F. Aguilera-Granja, 1 J. Ferrer, 2 and A. Vega 1 

1 Departamento de Fisica Teorica, Atomica y Optica. Universidad de Valladolid, E-47011 Valladolid, Spain 
^Departamento de Fisica. Universidad de Oviedo, Spain 
(Dated: February 1, 2008) 

We present a thorough theoretical assessment of the stability of non-collinear spin arrangements 
in small palladium clusters. We generally find that ferromagnetic order is always preferred, but that 
antiferromagnetic and non-collinear configurations of different sorts exist and compete for the first 
excited isomers. We also show that the ground state is insensitive to the choice of atomic configura- 
tion for the pseudopotential used and to the approximation taken for the exchange and correlation 
potential. Moreover, the existence and relative stability of the different excited configurations also 
depends weakly on the approximations employed. These results provide strong evidence on the 
transferability of pseudopotential and exchange and correlation functionals for palladium clusters 
as opposed to the situation found for the bulk phases of palladium. 

PACS numbers: 73.22.-f. 75.75.+a 



INTRODUCTION 

The magnetic properties of free-standing atomic clus- 
ters of 3d TM elements have been intensively scrutinized 
during the last two decades. Two different but related 
phenomena have specifically been discussed and essen- 
tially unravelled. The first is the modification of local 
magnetic moments as compared with the values found 
in bulk materials. The second is the competition be- 
tween the possible ferromagnetic, antiferromagnetic and 
non-collinear arrangements of the local spins, as well as 
its interplay with the geometry of the nanostructure. In 
the case of ferromagnetic elements like Fe, Co and Ni, 
the increase of the average cluster magnetic moment can 
be easily explained in terms of the reduced atomic co- 
ordination in the low-dimensional regime, with oscilla- 
tions associated to structural (symmetry) changes, [lj 
The case of antiferromagnets like Cr and Mn is much 
more complex. Atoms of these elements may display 
large magnetic moments, since they have a large num- 
ber of d-holes susceptible to be polarized. On the other 
hand, clusters of these atoms may display tiny average 
magnetizations due to the tendency of their atomic mo- 
ments to align in antiparallel directions. The structure 
plays also a fundamental role in the magnetic behavior 
of these clusters, since it may originate magnetic frus- 
tration. A conventional example of magnetic frustration 
in a classical spin system appears when atoms positions 
form triangular motifs. The studies of these classical sys- 
tems show that magnetic frustration frequently leads to 
non-collinear configurations of the local spin moments. 
The latest theoretical studies reported in the literature 
show that non-collinear arrangements of quantum spins 
also appear as the ground or as some of the first isomers 
of clusters of 3d atoms, including not only Cr and Mn, 
but also Fe, Co and Ni. 0, |, 1, [B[ @ 

All materials made of 4d TM elements are paramag- 
nets, in contrast to those of the 3d row. A natural ques- 



tion thus arises of whether small clusters of Ad elements 
may show low-lying magnetic states of collinear or even 
non-colinear nature. Bulk palladium, being a paramag- 
net in the brink of becoming a ferromagnet, presents one 
of the most intriguing and controversial magnetic behav- 
iors in nature. [7( It is therefore not surprising that the 
very few experimental and theoretical studies published 
so far try to clarify whether Pd clusters of given sizes are 
magnetic or not, and what is the order of magnitude of 
their average magnetic moment. From the experimental 
side, most of the reports agree that only very small clus- 
ters have a net magnetic moment SSIISJp, with the 
exception of Shinohara and coworkers. [l2| who found 
noticeable magnetic moments at the surface of Pd par- 
ticles as big as 79 A. From the theoretical side, there 
is also consensus that very small Pd clusters are indeed 
magnetic. [H [3, [U, OH, Gj, EH Futschek et al. [H have 
studied recently small Pd clusters using Density Func- 
tional Theory (DFT) in the collinear framework, within 
a fixed-moment mode. They have found that multiple 
spin isomers exist for each cluster size with very small 
energy differences. Interestingly, some of these compet- 
ing isomers present ferromagnetic order, while others dis- 
play antiferromagnetic alignments, with possible frustra- 
tion. Although Pd has tendency to ferromagentic order, 
this fact strongly points out to the possible existence of 
non-collinear magnetic structures, as a mechanism to re- 
lease the frustration and competition between the differ- 
ent magnetic solutions. 

We report in this article a thorough Ab initio study of 
the magnetic behavior of small palladium clusters Pdjv, 
with N ranging from 3 to 7. We have performed a si- 
multaneous optimization of the geometric and magnetic 
degrees of freedom fully allowing for non-collinear spin 
arrangements. A debate exists currently on the accu- 
racy of the Local Density Approximation (LDA) ver- 
sus the Generalized Gradient Approximation (GGA) [2l| 
for the determination of the magnetic behavior of low- 
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FIG. 1: Illustration of the ground state structures of the dif- 
ferent clusters here studied and average interatomic distances 
(in A) within GGA1. 
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The present 

article shows that both approximations provide the same 
results for the ground state of free-standing Pd atomic 
clusters. Moreover, even the existence and energy order- 
ing of the first excited states depend much more weakly 
on the choice of correlation functional and pseudopoten- 
tial in than the bulk fee material. 



THEORETICAL DETAILS 

We have performed our calculations using the code 
SIESTA. (Hi SIESTA is a DFT method that employs 
linear combination of pseudoatomic orbitals as basis 
set. The electronic core is replaced by a nonlocal 
norm-conserving Troullier-Martins[25| pseudopotential 
that may include nonlinear core correction terms. The 
code allows to perform, together with the electronic cal- 
culation, structural optimization using a variety of al- 
gorithms. It also allows to simulate non-collinear spin 
arrangements both in the LDA and in the GGA approx- 
imations. [26] 

In the present calculation, we have also used a variety 
of pseudopotentials to test their effect on free-standing 
clusters and their corresponding transferability. We have 
generated three different pseudopotentials using LDA. 
The first (LDAl) was built with the electronic config- 
urations 5s 1 , 5p° and 4d 9 , and core-corrections matching 
radius r c — 2.00 a.u.; the second (LDA2) was identical to 
LDAl, but with r c = 1.2 a.u.; the third had a closed-shell 
atomic configuration (5s°, 5p° and Ad 10 ) and r c = 1.2 a.u. 
We have also generated two GGA pseudopotentials with 
electronic configuration 5s 1 , 5p° and Ad 9 , and r c = 2.0 
or 1.2 a.u. (GGA1 and GGA2, respectively). In all five 
cases, the cutoff radii of the s, p and d orbitals were 
taken at 2.30, 2.46 and 1.67 a.u., respectively. We have 
described valence states by a double-^ polarized basis set 
(e.g.: two different radial functions for s and d orbitals 
and a single one for p orbitals). We have taken an energy 
cutoff of 150 Ry to define the real space grid for numer- 
ical integrations, but we checked that higher cutoffs did 
not alter the results. We have carried out the structural 
optimization using a conjugate gradient algorithm, where 
we have set the tolerance for the forces at 0.003 eV/A, 
with eventual double-checks using 0.001 eV/A. 



RESULTS 

We have found that the five pseudopotentials provide 
similar results when applied to an isolated palladium 
atom, being the eigenvalues of the ground state and 
different excited states slightly better reproduced with 
LDAl and GGAl (both had r c = 2.00 a.u.). However, 
we have observed that they give rise to different magnetic 
behaviors when applied to the bulk fee material. We 
first remind that palladium has a paramagnetic ground 
state with a lattice constant of 3.89 A. Moreover, the lat- 
est simulations show that LDA predicts a paramagnetic 
ground state with a lattice constant of about 3.85-3.91 
in agreement with experiments, while GGA predicts the 
ground state to be ferromagnetic, with a much larger lat- 
tice constant of about 4.0 A and a magnetic moment of 
about 0.4 [Lb- The fact that LDA provides a slightly 
better description of the ground state of bulk hep has 
been stressed recently by Alexandre and co-workers [23| . 
Our LDAl pseudopotential gives a ferromagnetic ground 
state with M w 0.54^i B , while LDA2 and LDA3 pre- 
dict the ground state to be paramagnetic, as it should. 
Finally, both GGA pseudopotentials lead to a ferromag- 
netic ground state with M ss 0.48/j.b. On the other hand, 
all LDA approximations give a lattice constant equal to 
3.90 A. In contrast, both GGA pseudopotentials predict 
it to be equal to 4.01 A. We note here that the fact that 
all LDA pseudopotentials predict the same lattice con- 
stant for the ground state regardless of its magnetic na- 
ture also happens with fee Iron, where LDA predicts the 
same ground state energy and lattice constant for the 
paramagnetic, antiferromagnetic and ferromagnetic Low 
Spin states 2a, 27 1. These results highlight the lack of 
transferability of the different pseudopotential and the 
sensitivity with the exchange and correlation functional 
used to describe the ground state of bulk palladium. 

It is therefore of great importance to assess whether the 
same difficulties apply to the case of the atomic clusters 
considered here. As we shall discuss below, atomic pal- 
ladium clusters are much more insensitive to the choice 
of pseudopotential and exchange and correlation func- 
tional. We shall show that the ground state is the same 
for the different approximations used here, in stark con- 
trast to the case of the bulk material. Further, a large 
portion of the low-lying excited states are reproduced by 
all approximations. 

Notice that we have not kept fixed the magnetic mo- 
ment in our simulations of the Pdjy clusters, but rather 
have allowed it to vary freely during the non-collinear 
iterative selfconsistency process, in contrast to previous 
authors. Moreover, while we can not rule out that we 
may have missed low lying solutions, we have endeav- 
ored to minimize this risk by feeding a large variety of 
non-collinear seeds for each cluster. This effort has al- 
lowed us to find a rich and complex family of metastable 
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TABLE I: Bindig energy of the ferromagnetic clusters in 
meV/atom. 



N 


LDA1 


LDA3 


GGA1 


Ref.flSl 


Ref.fl9] 


3 


1.755 


1.326 


1.289 


1.203 


1.250 


4 


2.293 


1.942 


1.769 


1.628 


1.675 


5 


2.502 


2.168 


1.933 


1.766 


1.805 


6 


2.721 


2.401 


2.110 


1.919 


1.949 


7 


2.791 


2.452 


2.155 


1 .953 


1 .985 



TABLE II: Inter-atomic distances of the ferromagnetic clus- 
ters in Angstrom. 



N 


LDA1 


LDA3 


GGA1 


Ref.flS] 


Ref.fl9] 


3 


2.49 


2.49 


2.56 


2.52 


2.52 


4 


2.57 


2.59 


2.64 


2.61 


2.61 


5 


2.61 


2.63 


2.70 


2.65 


2.64 


6 


2.62 


2.64 


2.70 


2.66 


2.66 


7 


2.64 


2.66 


2.72 


2.68 


2.70 



solutions, that was absent in previous works. We finally 
note that we have repeated all calculations with the pseu- 
dopotentials LDAl, LDA3 and GGA1. 

We have found that all clusters, except Pdg, share the 
same collinear magnetic ground state, with a total spin of 
2 Hb, 
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The 



in agreement with previous authors 
geometry of the ground state and the average interatomic 
distance of the Pdjv clusters is displayed in Fig. 1, where 
we show that these range from 2.55 A in Pd3 to 2.71 
A for Pd7. We have written the binding energies of the 
different clusters in Table I. The table shows that GGA1 
gives slightly smaller values than LDAl and LDA3, as 
otherwise expected. Moreover, the binding energies pre- 
dicted by GGA1 are very similar to those obtained by Ku- 
mar, who also used the GGA (within an ultrasoft pseu- 
dopotentials, plane waves code) and by Futschek et al., 
who used the all-electron VASP code, but did not state 
the approximation employed. 

We should stress that all the tested pseudopotentials 
provide the same ground state, in stark contrast to the 
situation that arose for the bulk material. Moreover, 
we have found very similar inter-atomic distances for all 
Pdjv clusters, as shown in Table II. The table shows that 
GGA1 predicts slightly longer distances than LDAl or 
LDA3, as otherwise. These distances also agree with 
those obtained by Kumar and Futschek within a range 
of 1 per cent. 

The Pd6 cluster displays a behavior different from the 
rest, and therefore we discuss it separately. Futschek 
and coworkers [l9[ found that Pde was also ferromag- 
netic in contrast to Kumar et al.[l8j, who predicted it 
to be paramagnetic. We have found that both states are 
nearly degenerate, with the paramagnetic solution being 
slightly more stable. Further, we have found that the 
relative energy of the two solutions show a strong depen- 



dence with inter-atomic distance. Fig. (2) shows that 
there is a level crossing when the bonds are elongated by 
just 1 per cent from the equilibrium distance. The close 
proximity between the equilibrium distance and the level 
crossing distance explains the sensitivity of these results 
with respect to slight differences in the simulations and 
the discrepancies between Kumar and Futscheck [l|| . 
Aditionally, we have been unable to find non-collinear or 
antifcrromagnetic solutions for this cluster. 

In contrast, and independently of the pseudopotential 
or approximation used, the rest of the clusters show a 
rich variety of antiferromagnetic and non-collinear solu- 
tions. Most of these solutions, though not all, exist for 
all LDAl, LDA3 and GGA1. We have also found that, 
whenever they exist, the relative order of the different 
solutions is maintained, and the size of the atomic mo- 
ments is very similar. These facts strengthen our belief 
that Pd atomic clusters are much more insensitive to the 
pseudopotential and approximation employed than bulk 
Pd. It is also reassuring that most of the collinear so- 
lutions have been identified in previous calculations [l9j| 
(e.g.: AF1 for Pd 4 and Pd 5 and AF2 for Pd 7 ). 

The non-collinear solutions found can be classified 
into those that release antifcrromagnetic frustration and 
therefore have lower excitation energy than the AF so- 
lution (NCI and NC2 in Pd 4 , shown in Table III and 
Fig. 3), and radial or quasi-radial solutions, that resem- 
ble the hedgehogs found in low dimensional theories of 
classical or quantum antiferromagnets 28] . Hedgehogs in 
these theories do not release frustration but rather are 
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FIG. 2: Energy difference of the Ferromagnetic state and the 
Paramagnetic one for the Pd6 clusters as a function of the 
average inter-atomic distance. We have taken the energy of 
the paramagnetic ground state as the zero of energy. Dis- 
tances are measured in units of the equilibrium distance of 
the cluster. The vertical dashed line indicates the theoretical 
interatomic distance of bulk palladium in the GGA1 approx- 
imation. 
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TABLE III: Different solutions obtained for the Pdjv clusters. We provide the absolute values of the atomic magnetic moments, 
the total magnetic moment in the cluster p, (both in units of fx B ) and the excitation energy per atom (in meV). For N=5 and 
7, the first two values of the atomic moments correspond to the axial sites, whereas the last ones correspond to the planar sites. 
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(0.67x3) 
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(0.67x3) 
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(0, 0.30, -0.30) 
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Radial 














(0.18x3) 





28 


Para. 


(0x3) 





68 


(0x3) 





34 


(0x3) 





75 


N=4 


Ferro. 


(0.50x4) 


2 





(0.50x4) 


2 





(0.50x4) 


2 





NCI 














(0.29,0.29,0.29,0.29) 





12 


NC2 


(0.35,0.24,0.24,0.35) 


0.25 
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(0.25,0.28,0.28,0.25) 


0.03 


10 








AF1 


(0.32, 0.32, -0.32, -0.32) 





26 


(0,23, 0.23, -0.23, -0.23) 





30 


(0.29, 0.29, -0.29, -0.29) 





25 


AF2 


(0.41, 0, -0.41, 0) 





40 


(0,31, 0, -0.31, 0) 





31 


(0.38, 0, -0.38, 0) 





36 


Para. 


(0x4) 





86 


(0x4) 





59 


(0x4) 





78 


N=5 


Ferro. 


(0.43,0.43,0.38x3) 
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(0.40x5) 
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AF1 


(0, 0, 0.43, - 0.43, 0) 
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(0, 0, 0.33, - 0.33, 0) 
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(0, 0, 0.39, - 0.39, 0) 
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AF2 


(0, 0, 0.48, -0.24, -0.24) 
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(0, 0, 0.44, -0.22, -0.22) 
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(0, 0, 0.29x3) 
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N=7 


Ferro. 


(0.19,0.19,0.32x5) 
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(0.21,0.21,0.31x5) 


2 





(0.20,0.20,0.32x5) 
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AF1 < 


^-0.36,0.36,-0.33,-0.22,0.22,0.32,0) 
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(0,0,-0.29,-0.20,0.20,0.29,0) 
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(0.27,0.27,0.18x5) 





22 


(0.24,0.24,0.12x5) 





14 


(0.24,0.24,0.17x5) 





20 


Para. 


(0x7) 





37 


(0x7) 





24 


(0x7) 





33 



excitations over the antifcrromagnctic ground state. We 
also find that these radial states have a higher energy 
that the antifcrromagnctic solution, and therefore do not 
release frustration. 

Notice that the antiferromagnetic and non-collinear so- 
lutions can be reached at temperatures of the order of 
room temperatures (25 meV). Therefore, any measure- 
ment of the magnetization performed at room temper- 
ature should find a thermal average of all those states, 
many of which have a tiny magnetic moment. It should 
not be surprising that such a measurement give a small 
net moment. 

We finally discuss the relationship between magnetism 
and equilibrium interatomic distances. We have found 
that these are essentially the same regardless of the mag- 
netic state for the largest clusters (N= 5-7), the smallest 
ones showing slight variations of less than 0.04 A, but 
only within the LDA solutions. We have additionally 



analyzed the relative stability of the different solutions 
as a function of the interatomic distance. To this aim, 
we plot in Fig. 4 the energy per atom of the low-lying 
excited states of the Pds cluster, relative to the ground 
state energy, as a function of an uniform volume expan- 
sion, obtained using GGA1. The figure shows that no 
crossover takes place, apart from the nearly-degenerate 
AF1 and AF2 solutions, that cross at an expansion of 
about 4%. Moreover, the relative energy differences are 
essentially preserved and the local magnetic moments 
kept constant, except for the AF2 and radial solutions, 
where they slightly change (by about 10%). 

CONCLUSIONS 

To summarize, we have studied the geometry and mag- 
netic properties of the ground state and lowest lying iso- 
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mers of small palladium clusters Pdjv, with N ranging 
from 3 to 7. Our results confirm that the ground state is 
indeed collinear or paramagnetic. We have found a rich 
variety of non-collinear low-lying isomers, some of which 
efficiently release frustration, while other (hedgehog-like 
solutions) do not. All these solutions should contribute to 
the room temperature magnetic behavior of the clusters, 
probably rendering small measured magnetic moments. 
We have finally found that simulations of atomic clusters 
are rather insensitive to the choice of the pseudopoten- 
tial and to the approximation used for the exchange and 
correlation potential. 

To provide a rough explanation of this contrasting be- 
havior, we first note that GGA and LDA provide results 
that are quantitatively very similar in most cases, includ- 
ing bulk palladium. GGA provides results that are qual- 
itatively different from LDA for bulk palladium because 
this material sits right at the edge of the paramagnetic- 
ferromagnetic transition, actually tilted slightly towards 
the paramagnetic side. Therefore, bulk palladium falls 
slightly short of fulfilling Stoner's criterion. GGA and 
LDA are mean field approximations. Hence, they tend 
to overstimate magnetic properties. This error is partly 
compensated by LDA, that predicts usually lower mag- 
netic moments than GGA. Consequently, LDA predicts 
the correct ground state for bulk palladium, while GGA 
provides a ground state that is slightly tilted towards the 
magnetic side of the phase diagram. Hence, LDA and 
GGA predict quantitatively similar, but qualitatively dif- 
ferent results for the magnetic properties of this material. 
Atomic Pd clusters are strongly magnetic, with the ex- 
ception of Pd 6 . Therefore, all these approximations pre- 
dict quantitative and qualitative similar results for those 
clusters. 
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FIG. 4: Excitation energy per atom of the magnetic solutions 
of the Pds cluster as a function of the average interatomic 
distance, using GGA1. 



ACKNOWLEDGMENTS 

This work was supported by the Spanish Ministerio de 
Educacion y Ciencia (grants MAT2005-03415, BFM2003- 
03156 and SAB2004-0129) in conjuction with the Euro- 
pean Regional Development Found, INTAS (Project 03- 
51-4778) and the Mexican (PROMEC-SEP-CA-230 and 
CONACyT-SNI). We wish to acknowledge useful conver- 
sations with J. M. Montejano-Carrizales, L. Fernandcz- 
Seivane, V. M. Garcia-Suarez and F. Yndurain. 




FIG. 3: Illustration of the non-colinear magnetic solutions 
for Pd 4 NC2 (LDA1), NC2 (LDA3) and NCI (GGA1). The 
arrows are proportional to the size of the atomic moments. 
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